Electron structure of the Falicov-Kimball model with a magnetic field 
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The two-dimensional Falicov-Kimball model in the presence of a perpendicular magnetic field 
is investigated by the dynamical mean-field theory. Within the model the interplay between elec- 
tron correlations and the fine electron structure due to the magnetic field is essentially emerged. 
Without electron correlations the magnetic field induces the electron structure to the so-called Hof- 
stadter butterfly. It is found that when electron correlations drives the metal-insulator transition, 
they simultaneously smear out the fine structure of the Hofstadter butterfly. In a long-range or- 
dered phase, the electron correlation induced gap preserves the fine structure, but it separates the 
Hofstadter butterfly into two wings. 

PACS numbers: 71.27. +a, 71.10.Fd, 71.70.Di, 67.85.-d 



I. INTRODUCTION 



The problem of electrons moving under an external 
magnetic field has attracted a lot of attention since the 
beginning of quantum mechanics. Two dimensional elec- 
tron gas under a perpendicular magnetic field creates the 
quantum Hall effect.— It comprehensively deals with 
the interplay of electron correlations and magnetic field. 
The integer quantum Hall effect is due to the quantiza- 
tion of the energy levels of free electrons under a magnetic 
fieldjii^ while the fractional quantum Hall effect is essen- 
tially due to the electron interaction under a magnetic 
fieldi^ The picture totally becomes complexity when 
electrons move additionally on a lattice or under a pe- 
riodic potential. The simultaneous presence of magnetic 
field and lattice potential gives the spectra of two di- 
mensional noninteracting electrons a fine structure of the 
famous Hofstadter butterfly^ The Hofstadter butterfly 
displays a recursive structure over rational gauge field 
and a Cantor set at any irrational gauge field. The Hall 
conductance of the noninteracting Bloch electrons is still 
quantized with an integer number when the Fermi en- 
ergy lies within a gap of the Hofstadter butterfly^ When 
electron correlations are included, the effect of simultane- 
ous presence of magnetic field and electron interaction on 
the lattice band energy remains interesting. Without the 
external field electron correlations can induce different 
phenomena, for instance, the metal-insulator transition, 
long-range ordered phases. When electron correlations 
are absent, the magnetic field creates the Hofstadter but- 
terfly of the electron structure. A simultaneous presence 
of electron correlations and magnetic field would induce 
a complexity of the electron structure. The exact diago- 
nalization of finite lattices shows that the electron corre- 
lations smear out the Hofstadter butterfly— However, it 
is not easy to distinguish the fine electron structure due 
to the presence of magnetic field from the discrete energy 
levels due to the finite-size effect of the exact diagonal- 
ization. The Hartree-Fock mean field calculations reveal 
the electron structure with the Hofstadter butterfly and 
an additional correlation-induced gap.SiE It gives rise to 
an interest in the study of the interplay between electron 



correlations and magnetic field in a two-dimensional lat- 
tice beyond the Hartree-Fock approximation and in the 
thermodynamic limit. In experimental aspect, with the 
rapid development of ultracold technique some funda- 
mental models of many quantum particles can be real- 
ized by loading ultracold particles into optical lattices 
(see, for example, the review in Ref. id 1 ). In particular, 
by using the technique of laser-assisted tunneling 11 or of 
lattice rotating^ artificial gauge fields can be created in 
optical lattices. As a results it is possible to realize the 
effect of magnetic field on the Bloch electrons by load- 
ing ultracold particles into optical lattices with artificial 
gauge field. Indeed, recently several models of optical 
lattices were proposed to study the Hofstadter butterfly 
of ultracold particles Jiii3iiiii^ 

In the present paper we theoretically study the effect 
of electron correlations on the Hofstadter butterfly of the 
electron structure under a magnetic field. The electron 
correlations are modelled by the Coulomb interaction of 
the Falicov-Kimball model (FKM)J£ It is a local repul- 
sive interaction of mobile electrons and massive localized 
particles. The FKM was originally introduced to describe 
a metal-insulator transition in transition-metal oxides. It 
can be viewed as a simplified Hubbard model where elec- 
trons with down spin are frozen and do not hop. The 
FKM was also used as a starting point to investigate dif- 
ferent electron correlation phenomena, for instance the 
mixed valence 1 ^ or the electronic ferroelectricity i 18 > 19 The 
FKM can also be incorporated into different models to 
study various aspects of electron correlations such as the 
charge ordered phase in manganese comp ound a 20 ' 2 1 ' 2 2 or 
electron localization. 23 ! 24 Much progress has been made 
on solving the FKM in both exact and approximation 
ways, where all properties of the conduction electrons 
are well known . 25 ' 26 ! 27 ' 28 In the homogeneous phase the 
FKM displays a metal-insulator transition. When the 
Coulomb interaction is strong, it prevents the mobility 
of itinerant electrons by forming the Mott-Hubbard gap. 
At low temperature a charge ordering occurs. For half 
filling the charge ordering gap opens at the Fermi level 
and drives the system into an insulating phase. One 
may expect that the electron-correlation induced gaps 
of the Mott-Hubbard type and of long-range charge or- 
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dering may have different effects on the fine structure of 
the Hofstadter butterfly. A realization of the FKM was 
also proposed as an optical lattice of a mixture of light 
fcrmionic atoms (e.g., 6 Li) and heavy fermionic atoms 
(e.g., 40 K)i 29 i 30 When the optical lattice modelling the 
FKM is established, it is also possible to create an ar- 
tificial magnetic field. In the present paper the two- 
dimensional square lattice with a perpendicular mag- 
netic field is considered. The dynamical mean-field the- 
ory (DMFT) is employed to calculate the electron struc- 
ture of the considered model. The DMFT is widely and 
successfully applied to study strongly correlated electron 
systemsi 31 ' 32 It gives the exact solutions in infinite di- 
mensions. However, for two-dimensional systems the 
DMFT is just an approximation. It neglects nonlocal 
correlations. However, the applications of the DMFT to 
FKM show that the approximation is still accurate in two 
dimensions. 33 ! 34 We find that the electron correlation ef- 
fect on the Hofstadter butterfly depends on the nature 
of the correlated phase when the magnetic field is ab- 
sent. When a long-range order is absent, electron correla- 
tions only induce the metal-insulator transition and they 
smear out the fine structure of the Hofstadter butterfly. 
In a long-range ordered phase, the electron correlation 
induced gap preserves the fine structure, but separates 
the Hofstadter butterfly into two wings. 

The plan of the present paper is as follows. In Sec. II 
we describe the FKM with a perpendicular magnetic field 
on a square lattice. We also present the DMFT for calcu- 
lating the Green function in this section. In Sec. Ill the 
numerical results are presented. Finally, the conclusion 
and remarks are presented in Sec. IV. 
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t ij = t exp ( i — J A • dl J , 
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where 0o = hc/e, and A is the vector potential. For a 
constant magnetic field perpendicular to the square lat- 
tice, the Landau gauge can be chosen for the vector po- 
tential A = (0, Bx,0), where B is the magnetic field 
strength. With this Landau gauge the hopping integral 
in the x direction is just t, while in the y direction it 
acquires additional phase factor t exp(±i2naxi), where 
a = Ba 2 /(j>Q, a is the lattice constant, and Xi is the lat- 
tice site position in the x-axis. In the following we will set 
a = 1. Parameter a is just the magnetic flux per unit cell 
in the units of the flux quantum 4>q. It is clear that the 
Hamiltonian is invariant with the translation a — > a + m, 
where to is a integer. Therefore it is only necessary to 
consider < a < 1. 

We will also only consider the rational magnetic field, 
i.e., a — p/q, where p, q are two coprime integers. The 
translation operator that moves q lattice spacing in the x 
direction leaves the Hamiltonian unchanged. We divide 
the lattice into q penetrating sublattice in x direction, i.c, 
each lattice site can be indexed by a number n, and its 
coordinates x, y, where n — mod(i? x , q), (1 < n < q). We 
take the Fourier transformation from the direct lattice to 
the reciprocal lattice 
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II. THE FALICOV-KIMBALL MODEL WITH A 
PERPENDICULAR MAGNETIC FIELD AND 
THE DYNAMICAL MEAN-FIELD THEORY 



In this section we present the DMFT for the Falicov- 
Kimball model in the presence of a magnetic field. The 
model is described by the following Hamiltonian 

H = £ t ijC \ Cj - /i ]T c\ Cl +E f J2 flfi 

+uj24ciflfi, (i) 



where cj (ej), fj (fi) are the creation (annihilation) op- 
erators for itinerant and localized electrons at site i, re- 
spectively, tij is the hopping integral of itinerant elec- 
trons between site i and j. U is the local interaction of 
itinerant and localized electrons, fi is the chemical po- 
tential for itinerant electrons and Ef is the energy level 
of localized electrons. We will consider only the half fill- 
ing case, where /i = — Ef = U/2. In the presence of a 



where N is the number of lattice sites. The wave vectors 
k x , k y are restricted to the reduced Brillouin zone 
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The hopping part of Hamiltonian in Eq. |T]) can be rewrit- 
ten as 



H t = J2xtE(k)X k , 



(3) 
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magnetic field the hopping integral acquires the Peierls G(k, w) = ((A\|A^}}w = [uj + p — E(k) — S(w)] , (5) 



with e,^ = cos(k y + (n — l)2irp/q). 

We apply the DMFT to the calculation of the Green 
function of itinerant electrons in the matrix form 

t\ 



where is the self energy. Within the DMFT the self 
energy is independent on momentum. Moreover, it is also 
diagonal, i.e., Ti nm (u>) = d nm Yi n (u>). This formulation is 
similar to the DMFT applications for the antiferromag- 
netic or checkerboard charge ordered phases^ Basically, 
the DMFT is exact in infinite dimensions. However, its 
application for two-dimensional systems is just approx- 
imation. The approximation neglects nonlocal correla- 
tions which exist as the momentum dependence and off- 
diagonal elements of the self energy. The DMFT calcu- 
lations for two-dimensional FKM without the magnetic 
field show that the approximation still accurate for the 
electron dynamicsj 33 i 34 

The self energy S„(w) is self consistently determined 
from the dynamics of a single interaction site embedded 
in an effective mean-field medium. For the FKM the 
effective single-site problem can be solved exactly*^ We 
obtain the Green function of the single site problem^ 
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where G n (u)) is the Weiss field for sublattice n. The 
weight factors W n0 and W„i can be calculated from the 
Weiss field 



W nl = f(E n ), 

W n0 = i-w nU 



(7) 
(8) 



where f(u>) = l/(exp(w/T) + 1) is the Fermi-Dirac dis- 
tribution function, and 
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FIG. 1: (Color online) The density plot of the DOS of itin- 
erant electrons of the high temperature phase (T = 1) for 
various values of U. The magnetic field parameter a = p/q 
varies with p — 0, 1, . . . , q, and q — 40. The color scheme 
corresponds to the rainbow color scheme, i.e. the red (violet) 
color corresponds to the largest (smallest) value of the DOS. 
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The Weiss field Green function Q n {w) also satisfies the 
Dyson equation of the effective single site problem, i.e., 



(10) 



The self-consistent condition requires that the Green 
function obtained from the effective single-site problem 
must coincide with the local Green function, i.e., 
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With this self-consistent condition the system of equa- 
tions for the self energy is closed. We can solve the sys- 
tem of equations by iterations as usuali^ 



III. NUMERICAL RESULTS 

In this section we present the numerical results ob- 
tained by solving the DMFT equations by iterations. We 
take t = 1 as the energy unit. The magnetic parameter 
a = p/q varies with p = 0,1,..., q. We take q = 40. 
When p and q are not coprime, we can reduce them to 



coprime integers. This also reduces the matrix dimension 
of the Green function, and saves the computation time. 
First we consider the high temperature phase, where the 
metal-insulator transition of the Mott-Hubbard type oc- 
curs. We emphasize that in this phase temperature doest 
not affect on the electron structure. The phase is just 
the homogeneous solution of the DMFT equations and 
it is stable at high temperature. The electron structure 
can be imaged by using the density plotting of the den- 
sity of states (DOS) of the itinerant electrons. In Fig. [T] 
we plot the image of the DOS for various values of U. 
It shows that the electron structure mimics the Hofs- 
tadter butterfly when the electron correlations are in- 
cluded. For weak interactions the fine structure of the 
Hofstadter butterfly still survives. However, the elec- 
tron correlations already smear out it. Some fine gaps in 
the structure of the Hofstadter butterfly are closed. As 
the value of U increases, the smearing becomes stronger. 
For strong interactions all fine gaps are closed. How- 
ever, a middle rough gap opens for U > U c r* 4£. This 
gap is essentially the Mott-Hubbard gap, which opens 
in the insulating phase. Without the magnetic field the 
metal-insulator transition occurs at U c . In the presence 
of the magnetic field the metal-insulator transition still 
occurs, but the lower and upper bands mimic the Hof- 
stadter butterfly. The electron structure is symmetry 
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FIG. 2: The DOS of itinerant electrons of the high tempera- 
ture phase (T = 1) for various values of a and U. 

in respect to lines u> = and a — 1/2 like the nonin- 
teraction case, i.e., p(uj,a) = p(—uj,a) = p(uj, 1 — a), 
where p(u>, a) = — lmG n (uj)/qir is the DOS of itiner- 
ant electrons. 

In Fig. [5] we plot the DOS of itinerant electrons for 
various values of a and U. When the electron correla- 
tions are absent (U = 0) the number of bands is just 
g£ One can observe that the q bands can be grouped 
into subgroups of bands which are separated by moder- 
ate gaps. Within a band subgroup the bands are also 
separated by fine gaps. For example, when a = 3/8, 
there are 3 subgroups of bands, one at w = 0, and two 
others around lo — ±2. When the interaction is included, 
the band number is reduced by closing the gaps. As 
the interaction increases, first the fine gaps within the 
band subgroups are closed, and then the gaps between 
the band subgroups are closed. In the insulating phase 
all the gaps in the Hofstadter butterfly are closed, but 



U-l U-4 




FIG. 3: (Color online) The density plot of the DOS of itiner- 
ant electrons of the charge ordered phase at low temperature 
(T = 0.01) for various values of U. The magnetic field param- 
eter a — p/q varies with p = 0, 1, . . . , q, and q — 40. The color 
scheme corresponds to the rainbow color scheme, i.e. the red 
(violet) color corresponds to the largest (smallest) value of 
the DOS. 



FIG. 4: (Color online) The DOS of itinerant electrons of the 
checkerboard charge ordered phase at low temperature (T = 
0.01) for various values of a and U. The red and green color 
lines correspond to the DOS of two penetrating sublattices of 
the checkerboard charge ordered phase. 



the Mott-Hubbard gap opens. The metal-insulator tran- 
sition occurs at the same value of U for all values of the 
magnetic field. Despite of closing of the fine gaps in the 
insulating phase, the intensity of the DOS of itinerant 
electrons still shows a mimic Hofstadter butterfly with 
smearing out fine features, as shown in Fig.[TJ 

Without the magnetic field the FKM also displays the 
checkerboard charge ordering at low temperature for any 
interaction U ^ 0. Wc study the possibility of the charge 
ordering when the magnetic field is present. In this case 
we additionally divide the lattice into two penetrating 
sublattices in the y direction. If q is odd integer, we 
double the value of p and q that a — 2p/2q and the 
checkerboard charge ordering is commensurate with the 
magnetic structure. In Fig.[3]we plot the image of the to- 
tal DOS of itinerant electrons in the checkerboard charge 
ordered phase. It shows that the fine structure of the Hof- 
stadter butterfly still remains, however, it is separated 
by a middle gap. The middle gap is the charge ordering 
gap, which locks itinerant electrons into the checkerboard 
charge pattern. As the interaction increases the charge 
ordering gap also increases, and the band width of the 
lower and upper bands are slightly reduced. The preser- 
vation of the fine structure of the Hofstadter butterfly 
was also observed within the Hartrec-Fock mean field 
calculations. 8,9 However, the Hartree-Fock mean-field ap- 
proximation cannot find the smearing of the Hofstadter 
butterfly due to electron correlations when a long-range 
order is absent. 

In Fig. 2] we also plot the DOS of itinerant electrons 
of the two penetrating sublattices of the checkerboard 
charge ordered phase for various values of U and a. It 
shows that the magnetic field does not affect on the 
charge ordering gap. The gap sorely depends on the in- 
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teraction as in the case of absence of the magnetic field. 
When q is even, the number of bands is also q, like the 
noninteraction case. However, unlike the noninteraction 
case, the subgroup of bands at the Fermi level is sepa- 
rated by the charge ordering gap. When q is odd, the 
band at the Fermi level is split into two bands which are 
also separated by the charge ordering gap, as shown in 
Fig. Q] (the case a = 2/5). The gaps between subgroups 
of bands are slightly reduced as the interaction increases. 

IV. CONCLUSION 

In the present paper we study the effect of electron cor- 
relations on the Hofstadter butterfly which is the electron 
structure of the two-dimensional Bloch electrons under a 
perpendicular magnetic field. By employing the DMFT 
we calculate the Green function of itinerant electrons 
in the case of rational magnetic field. Electron correla- 
tions exhibit different effects on the Hofstadter butterfly 
depending on the nature of the correlated phase when 
the magnetic field is absent. In the absence of a long- 
range order, electron correlations also induce the metal- 
insulator transition when the magnetic field is present. 
However, the electron correlations smear out the fine 
structure of the Hofstadter butterfly. The number of 
bands is reduced as the interaction increases. In the in- 
sulating phase, all fine gaps of the Hofstadter butterfly 



are closed, but the Mott-Hubbard gap opens at the Fermi 
level. In a long-range ordered phase, the electron correla- 
tion induced gap, such as the checkerboard charge order- 
ing gap in the FKM, preserves the fine structure of the 
Hofstadter butterfly. However, the Hofstadter butterfly 
is separated into two wings by the long-range ordering 
gap. 

In the present paper we have only considered the ra- 
tional magnetic field. In the noninteraction case an ir- 
rational magnetic field induces the Hofstadter butterfly 
in the form of a Cantor set. The magnetic structure is 
incommensurate with the lattice structure as well as with 
the checkerboard charge ordering pattern. The effect of 
electron correlations on the such Hofstadter butterfly re- 
mains open, and we leave it for further study. 
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